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Abstract We consider a stochastic, continuous state and time opinion model 
where each agent’s opinion locally interacts with other agents’ opinions in the 
system, and there is also exogenous randomness. The interaction tends to create 
clusters of common opinion. By using linear stability analysis of the associated 
nonlinear Fokker-Planck equation that governs the empirical density of opinions 
in the limit of infinitely many agents, we can estimate the number of clusters, the 
time to cluster formation and the critical strength of randomness so as to have 
cluster formation. We also discuss the cluster dynamics after their formation, the 
width and the effective diffusivity of the clusters. Finally, the long term behavior 
of clusters is explored numerically. Extensive numerical simulations confirm our 
analytical findings. 
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1 Introduction 

Opinion dynamics models have attracted a lot of attention and there are many 
analytical and numerical studies that consider different models arising from many 
different fields. In much of the literature, an opinion dynamics model is a system 
with a large number of opinion variables, xi(t), i = 1 ,... ,N, taking values in R". 
The time evolution of the opinion variables is governed by an attractive interaction 
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between any two opinion variables, often taken to be a nonnegative function of 
the Euclidean distance of the two opinion variables and may also be time depen¬ 
dent. The most interesting feature of such a model is that opinions only interact 
locally and the influence function is compactly supported, interpreted as bounded 
confidence. In this case, it is of interest to know whether the system will exhibit 
consensus convergence , which means that all the opinion variables converge to the 
same point as time tends to infinity. Except for some specific consensus models, 
a broad sufficient condition to have consensus convergence for a general class of 
models is not known. However, several studies have shown that for a variety of 
different types of consensus building interactions, and without external forces or 
randomness, the opinions will converge to possibly several clusters. In this case, 
the distance between distinct clusters should be larger than the support of the 
influence region. But it is not known, in general, how to determine the number of 
clusters. 

A more realistic way to model opinion dynamics is to add external randomness 
to the system. In this case, the model becomes a system of N stochastic processes 
and usually the randomness in the model is independent from one opinion holder or 
agent to another. Many deterministic techniques can also be used in the stochastic 
case, but some methods, such as the use of master equations, are particularly useful 
in stochastic models. When the external noise is large in the stochastic models then 
the tendency to cluster is effectively eliminated as the system is dominated by the 
noise. This is a phenomenon seen elsewhere in statistical physics as well. The 
strength of the noise or randomness must be below a critical value in order for 
cluster formation to emerge and evolve. 

The literature in opinion dynamics is very extensive so we mention only a few 
papers that have guided our own work. Hegselmann and Krause m consider a 
discrete-time evolution model, in which the opinions in the next step are the aver¬ 
age of the current opinions within a specified range of the influence region. Pineda 
et al. [25] add noise to the Hegselmann-Krause model and determine the critical 
strength of the noise so as to have cluster formation, using a master equation ap¬ 
proach and linear stability analysis. The same method is also used in |23j |24| on 
the Deffuant-Weisbuch model 0-In a, the authors take the limit as the number 
of opinions goes to infinity and consider the distribution of the opinions (the Eule- 
rian approach), instead of tracking every single opinion in the Hegselmann-Krause 
model (the Lagrangian approach), and [21] further discuss the case with external 
forces. The long time behavior and a sufficient condition for consensus convergence 
of the Hegselmann-Krause model are considered in BUSIES] . The long time behav¬ 
ior of the Hegselmann-Krause model with a general influence function is discussed 
in 116] [22] . The Hegselmann-Krause model involving different types of agents is 
considered in m- Some recent development of the study of opinion dynamics are 
in mm- Other, related relevant works are BH PH EH] El BD El BH BH] B! - 

Our contributions in this paper are the following. We consider a stochastic 
opinion model where every opinion is influenced by an independent Brownian 
motion. By the mean field limit theory, the empirical probability measure of the 
opinions converges, as the size of the population goes to infinity, to a solution of 
a nonlinear Fokker-Planck equation. Using a linear stability analysis, we estimate 
the number of clusters, the time to cluster formation and the critical strength 
of the Brownian motions to have cluster formation. The linear stability analysis 
can be applied to both deterministic and stochastic models. We also discuss the 



Consensus Convergence with Stochastic Effects 


3 


behavior of the system after the initial cluster formation but before further clus¬ 
ter consolidation, where the centers of the clusters are expected to behave like 
independent Brownian motions. Finally, we consider the long time behavior of 
the system. Once clusters are formed, their centers behave like Brownian motions 
until further merging. After consensus convergence, where there is only one clus¬ 
ter, there is a small probability that all the opinions inside the limit cluster will 
spread out and the system will become an independent agent evolution. Extensive 
numerical simulations are carried out to support our analysis and remarks about 
cluster formation and evolution. 

The paper is organized as follows. The interacting agent model is presented 
in section [2] The mean field limit is presented briefly in section [3] The linearized 
stability analysis of the governing nonlinear Fokker-Planck equation is presented in 
section 0 ] when there is no external noise. The results of numerical simulations are 
also presented in this section. In section [5] we extend the analysis of the previous 
section to the stochastic case when there are external noise influences. We also 
present the results of numerical simulations in the stochastic case. In section [6] we 
comment briefly about the long time behavior in the stochastic case when there 
is clustering. We end with a brief summary and conclusions in section [ 7 ] 


2 The interacting agent model 

The opinion model we consider in this paper is (see (221 Ecp (1.2a)]): 

1 N 

dxi = Y2 a ij( x i - x j) dt + a dW l (t), (1) 

j=i 

where Xi(t) is the agent V s opinion modeled as real valued process, where t is 
time and i = 1,..., N. The coefficients denote the strength of the interaction 
between and xj and they are a function of the distance between Xi and Xj\ 

CLij = <j)(\Xi- Xj\). (2) 

The interesting case is when the influence function <j> is non-negative and com¬ 
pactly supported. In other words, the interactions are attractive and the agent i 
affects only the other agents that have similar opinions. Here we assume that (j> is 
compactly supported in [0,1?o] 

Hr) = 00 (^) (3) 

where supp(()>o) = [0,1]. 

The W l (t ), i = 1,..., TV are independent standard Brownian motions that 
model the uncertainties of the agents’ opinions, and a is a non-negative constant 
quantifying the strength of the uncertainties. If a s= 0, then there is no randomness 
in this model and 0 is a deterministic system, while if a > 0, the system becomes 
stochastic. 

For the purposes of the analysis below, we consider the model 0 on the torus 
[0, L\ instead of the real line R. i.e. we consider the model in the bounded space 
[0, L] with periodic boundary conditions. The assumption of periodic boundary 
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conditions is mostly for simplifying the analysis. Although this assumption may 
not be appropriate in some applications, we found that the results obtained using 
it are numerically consistent with the same model in full space or in a finite interval 
with reflecting boundary conditions. The later two are in general more realistic for 
many applications. We note that the same periodicity assumption for the analysis 
of the opinion dynamics is also used in E3H2UBS]. 


3 The mean field limit 


At time t, we consider the empirical probability measure p N (t,dx) of the opinions 
of all the agents: 

1 N 

p N (t, dx) = — S Xj (t) (dx). (4) 

3=1 

Here S x (dx) is the Dirac measure with the point mass at x. The empirical proba¬ 
bility measure p N (t, dx ) is a measure valued stochastic process. We assume that as 
N —> oo, p N ( 0, dx) converges weakly, in probability to po(dx) which is a determin¬ 
istic measure with density po(x). By using the well known mean held asymptotic 
theory (see, for example, [T] [TQj [lj]), it can be shown that as N —> oo, p N {t,dx) 
converges weakly, in probability to p(t,dx ), for 0 < t < T < oo, a deterministic 
probability measure. Under suitable conditions the limit measure has a density 
p(t, x) which satisfies (in a weak sense) the nonlinear Fokker-Plank equation: 


dp 

dt 


(t,x) 


d_ 

dx 


I 


p(t,x- y)y<j>(\y\)dy 


, ,1 a 2 d 2 p, , 

+ ~ 2 w*)’ 


( 5 ) 


with given initial density po(x). In particular, if a:i(0),..., xjv(O) are sampled in¬ 
dependently and identically according to the uniform measure over [0, L], then the 
result holds true and the initial measure has constant density po(x) = 1/L. 

In this paper, we assume that N is large and view the mean held limit as 
the dehning problem. Therefore, we will analyze the overall behavior of the opin¬ 
ion dynamics, x% (t), - • • by analyzing instead the nonlinear Fokker-Planck 

equation <©■ 


4 Deterministic consensus convergence: er = 0 

We will follow a modulational instability approach to study the mean held limit 
when (cr = 0), also analyzed in [22] [16]. We look for conditions so as to have 
consensus convergence where all the opinions converge to a cluster as t —> oo. We 
also analyze the number of clusters if there is no consensus convergence and the 
time to cluster formation, that is, the onset of cluster formation. 


4.1 Linear stability analysis 

We hrst linearize the Fokker-Planck equation © with a = 0 by assuming that 
p(t,x) = po + pi(t, x) = 1/L + p\(t,x). Substituting p(t,x) = po + pi(t,x) into 
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([H]) and assuming that pi is a small perturbation of p so that the 0 (p\) term is 
negligible, we find that p\ ( t , x) satisfies: 


dpi 

dt 


(t,x) 


J Poy<l>(\y\)dypi(t,x) 

dx 

j Pi {t,x-y)y<j>{\y\)dyp 0 

J \) dy - 



( 6 ) 


The last equality in © holds because (p(\y\) is an even function and therefore 

JyH\y\) d y = 0 . 

By taking the Fourier transform in x, pi(t,k) = pi(t, x)e~ lkx dx, with the 
discrete set of frequencies k in 


K, = { 2nn/L , n £ Z}, 

we find from © that 

^-( t,k)= ip 0 k J e~ lkv y<t>(\y\)dy 
which gives the growth rates of the modes: 


Pi(t,k), 


7 k = Re 


ipok J e lfc!/ #(|y|)ch/j = p 0 k J sm(ky)y<j>(\y\)dy. 


( 7 ) 

( 8 ) 

( 9 ) 


We can see that for each k , \pi(t, k)\ = |pi(0, k)\ exp (7 ^t). By replacing 0 with 0o 
(see ©)> we can rewrite 7 ^ as 


Ik = poRo^(kR 0 ), 



4>o(s)s sin (qs)ds. 


( 10 ) 


The growth rate 7 ^ is maximal for k = ±fc ma x with fc m ax :— g ma x/-Ro, more exactly, 
for k equal to plus or minus the discrete frequency fc ma x in the set K. that maximizes 
tp(kRo), which is close to g ma x/-Ro- Here 


g ma x = argmax [-i/>(g)] = argmax 
q> 0 q> 0 



sin(gs)ds 


( 11 ) 


The optimal (continuous) frequency g ma x is positive and Hnite under general con¬ 
ditions since ~ 2 g 2 J 1 s 2 4>o(s)ds for 0 <g<l and V'(g) is bounded or decays 
to zero at infinity depending on the regularity of <f> q. 


4.2 Fluctuation theory 

By the central limit theorem, if we assume that the initial opinions *i(0),..., xjy(0) 
are sampled identically according to the uniform distribution over the domain 
[0, L\, then 


N 

Pi 


— ( N 
(t = 0 ,dx) := VN p N (t = 0 ,dx) - po{dx)^j = y/N I 77 ^<^(0 )(dx) 

V i =1 


dx 

T 
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converges in distribution as IV —► oo to the measure pi(t = 0, da;), whose frequency 
components, for k £ ZC\{0}, 


Pi{t = 0 ,k) 



h N 

I NJ2 S *i(0)( dx ) 



= lim VN 

N —»-oo 


( 1 N 

I ir e -fcj(°) 

\ 3= 1 



= lim VN 
N—too 



—ikxj (0) 


are independent and identically distributed with complex circular Gaussian ran¬ 
dom variables with mean zero and variance 1: 

E \pi(t = 0, k)] = 0, E [pi (t = 0, k)pi (t = 0, A/)] =4fc', k,k'eK C\{0}, 


pi(t = 0, —k) = pi(t = 0, k), while pi{t = 0, k = 0) = 0. 

For any T, the measure-valued process 

pi (t,dx) := \/iV (p N (t,dx) — po(dx)) , t£[0, T] (12) 

converges in distribution as JV -> 00 to a measure-valued process pi (t, dx) whose 
density pi (t, x) satisfies the deterministic PDE 

^-(t,x)t=p 0 J ^-(t,x-y)y<j>(\y\)dy (13) 

with the random initial condition described above. 

Consequently, combining with © and ©, at any time t, the frequency com¬ 
ponents pi(t,k), k £ AT\{0}, are independent complex circular Gaussian random 
variables, with mean zero and variance exp( 27 ^t): 


E ^pi(t,k)pi(t,k')^ = S kk > exp( 27 fc t), k,k'eJC\{ 0}, (14) 


pi(t, —k) = pi(t, k), while pi(t, k = 0) = 0. Therefore, 


E [p 1 (t,x)pi(t,x')] = E 


^2pi(t,k) 

. k 


gikx 

~zT 


k 


gikx' 

L 


= E E 


k 


gikx g— ikx' 

h (t, (*’ ~ k ) L 


J_ g 27 k t e ik(x-x') 

L 2 

k^O 


For large times, the spectrum of pi(t, x) becomes concentrated around the 
optimal wavenumber k m ax- We can expand 7 /. = 7max + 37 max(^ — Zc ma x) 2 for k 













Consensus Convergence with Stochastic Effects 


7 


around fc max and use a continuum approximation for the discrete sum: 

E [pi (t,x) Pl (t,x')\ e 2 ™V fc(x - x ') ~ ± f°° e^e^-^^dk 

k^O 2 ~°° ^ 


1 

27tL 


f 

—Re I e 27maxt [ 

L \ Jo 


7 lL 


27 fetg ik (. x ~x') c lk ~ —Re (| e 2(7m,x+i7^„(*:-km»x) 2 )t e !l;(x-i') 


dfc 




7tZ/ 

1 


= l r 


2'Ymax t 


cos(fc max (a; - a/)) 


\/ 717m a x 1 1 - 


dk 


_ (x — x ') 2 

e 4 l'l'maxl‘ 


(15) 


A typical realization of pi (t, x) is a modulation with the carrier spatial fre¬ 
quency fcmax and a slowly varying envelope with stationary Gaussian statistics, 
mean zero, and Gaussian covariance function. This is valid provided L 2 4|7" ax |t. 
If L 2 -C 4|7m ax |£, then the continuum approximation is not valid and we have 


E [p 1 (t,x)pi(t,x')\ m ~ ^2 e 27fc V fe(x * ] = J2^2 e 2lkt cos (k(x - x)) 


k^o 


L 2 


k> o 


- L 2 e 


27„ 


' cos (fc max (a; - x)) . 


(16) 


A typical realization of pi (t, x) is a modulation with the carrier spatial frequency 


Because 7max > 0, the linear system ( |13| ) is unstable and therefore the central 
limit theorem cannot be extended to arbitrarily large times. In fact the theorem 
is limited to times t such that pi(t,x)/VN is smaller than po = 1/L so that the 
linearization around po is valid. Therefore the time up to the onset of clustering 
is when the perturbation pi becomes of the same order as %/Apo, that is to say 
when E[pi(t c i u ,a;) 2 ] ~ AL~ 2 , which is approximately (up to terms smaller than 
In A): 

tclu - — In A ~ -—— J -v In A 

^7max Wmax) 


when A » 1. 

We note that the fact that a random initial distribution gives rise to a quasi- 
deterministic subsequent evolution by spectral gain selection occurs in many fields, 
for instance in fluid mechanics (hydrodynamic instabilities) or in nonlinear optics 
(beam hlamentation). 


4.3 Consensus convergence 


The linear stability analysis shows that the opinion dynamics, starting from a 
uniform distribution of agents, gives clustering with a mean distance between 
clusters equal to 2n /fc max (see (15) and |l6|)). Once clustering has occurred, two 
types of dynamical evolutions are possible: 


1. If 27r/fc max > Ro, then the clusters do not interact with each other because 
they are beyond the range of the influence function. Therefore, the situation is 
frozen and there is no consensus convergence. 
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2. If 27r/fc ma x < Ro, then the clusters interact with each other. There may be 
consensus convergence. However, consensus convergence is not guaranteed as 
clusters may merge by packets, and the centers of the new clusters may be sepa¬ 
rated by a distance larger than 27r/fc ma x, and then global consensus convergence 
does not happen. The number of mega-clusters formed by this dynamic is not 
easy to predict. 


If we neglect the rounding and consider fc max = Smax/Ho, which is possible if 
<?maxT/-Ro ^ 1 , then the criterion 2n/k ma .x > or < Rq does not depend on Ro, as it 
reads 27r/g max > or < 1, which depends only on the normalized influence function 
4>o by (10) and (111. 

These two dynamics can be observed in the examples of Figure 1.1 in [22] : 

1. If <j>{r) = l[ 01 ](r), then g max — 2.75 and the mean distance between clusters is 
about 2.3, that is beyond the range 1 of the influence function, and there is no 
consensus convergence. 

2. If <p(r) = 0.1 x 1[ Q ( r ) + l[i/ ^ i] ( r ), then q max ~ 9.1 and the mean distance 

between clusters is about 0.7, that is within the range of the influence function, 
and there is consensus convergence. 


These predictions are quantitatively in very good agreement with the numerical 
simulations (distance between clusters and so on). 

To summarize, the main result in the noiseless case a = 0 is as follows. In 
the regime N —> oo, there is no consensus convergence if g max < 27r. There may 
be consensus convergence if g ma x > 2-7T. Of course this stability analysis and the 
result that follow can be extended easily to the multi-dimensional case, and to 
other types of opinions or flocking dynamics. 


4.4 Numerical simulations 

We use the explicit Euler scheme to simulate the deterministic opinion dynamic 
(JT]) when <j = 0: 

1 N 

xl +1 =Xi- - x'jDixi -Xj)At, = c/>o{s/Ro). (17) 

3=1 

Although our analysis is on the torus [0,L], we still simulate © on the full 
space. The simulation results indicate, however, that the analysis under periodic 
conditions is still consistent with the numerics with different boundary condi¬ 
tions. As it is shown in [22] [16], if *i(0),..., *at( 0) are in the interval [0, L] . then 
x\ (t),..., xpf(t) £ [0, L\ for any t > 0. 

We test for the influence functions studied in |22| 1 16] : 

^o( s ) = 1 [0,l/V2]( s ) + 0- 1 X 1 (1 / n /2,1](«) 

00 (s) = l[o,i](s) 

00 (s) = 0.5 X l [0jl/v /2](s) + 1 (1 / n /2,1](«) 

00 ( s ) = 0- 1 X 1 [0,l//2]( S ) + 1 (l/v / 2,l] ( S ) 

0o( s ) = (1 ~ s ) 3 x l[o,i]( s ) 

0o( s ) = (1 — s ) 6 x 1 [ 0 , 1 ]( S ); 





(f>o{s) (t>o{s) 0o (s) 
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Fig. 1 The plots of the influence functions rf>o(s). 


and their plots are shown in Figure [l] 

We compute the key quantity g ma x by exploring all possible q in [0,100]: 


<7max = argmax 

R o q&lC,0<q<100 



sin (qs)ds 


We find that for the cases of and <j> g, g m ax are not unique and the non-uniciueness 
of <z ma x will greatly affect the results of the consensus convergence. The parameters 
we use for the simulation are At = 0.1, L = 10, Rq = 1 and N = 500. For each <f>o, 
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M s ) — 1 x l[o,l/\/2] ( s ) T 0-1 x l(l/-v/2,l] ( s ) 

A'' = 500, R 0 = 1, ?m»x = 3.77 


Inter-Cluster Distance = 1.67 



Time, t c i u = 21 


Fig . 2 Simulation for </>o(s) = 4>q(s). Left: 'ifj(q) evaluated at Roq G JC. Right: Simulations of 
( |1T[ ). The vertical dashed line is at t = t c i u . 


<A) (s) = 1 X l [ 0 ,i/^] W + lx 1 (l/x/2,1] ( s ) 



0o(s) = 1 X 1 |0 . 1/V 2] (s) + 1 X (s) 

IV = 500,-Ro = 1, 9max = 2.51 



Time, t du = 14.9 


Fig . 3 Simulation for <^>o( s ) — 0o( s )- Left: if>(q) evaluated at Roq G JC. Right: Simulations of 
The vertical dashed line is at t = t c i u . 


we also plot the function ip(s) = 2 q f Q c/>o(s)ssin(qs)ds ; the stars in the plots are 
the values of ip(s) at Roq € AC and the lines are the continuum approximation. 

From Figure[2j we can see that there is a unique g ma x = 3.7699. From our anal¬ 
ysis, we do not expect to have consensus convergence because g ma x = 3.7699 < 27r. 
The distance between clusters is 27ri?o/(j , m a x = 1.67, and therefore we should have 
roughly L/1.67 = 5.99 clusters, and indeed we have 6 clusters in our simulation. 
In addition, t c i u = In N/(2poRoi>(qina.x)) = 21 (the vertical blue dashed line) also 
correctly predicts the time to cluster formation. 

In Figure[3j if <j>o(s) = 1 [o.i]( s ), then ip(q) has a unique q ma x = 2.51 but it also 
has many suboptimal q where tp(q) is very close to ip(qma.x)- Note that q m ax < 2-7T 
means that there is no consensus convergence. The inter-cluster distance is 2.5, 
which is correct for the top and the bottom clusters. However, the central clusters 
are affected by the suboptimal q and therefore their inter-cluster distances are 
different. We also note that t c i u = lnN/(2poA7oV’( ( i , m a x)) = 14.9 (the vertical blue 
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0o( s ) — 0-5 X ^-[0,l/\/2] ( s ) + 1 x l(l/v/2,l] ( s ) 



</>o(s) = 0.5 X l [0>1/ ^ 2 ](s) + 1 x l {1/v ^ il] (s) 
N = 500, R 0 = 1, gmax = 53.4 
Inter-Cluster Distance = 0.118 



Tim e,t du = 11.5 


Fig . 4 Simulation for </>o(s) = ^o( s )- Left: 'ifj(q) evaluated at Roq E /C. Right: Simulations of 
( |1T[ ). The vertical dashed line is at t = t c i u . 


M s ) ~ 0-1 x l[o,i/x/2] ( s ) + 1 x l(i/v^,i] ( s ) 



M s ) — 0-1 X -*-[0,l/\/2] ( S ) + 1 X ^-(l/V2,l] ( S ) 
N = 500, R{) = 1,9 max = 53.4 


Inter-Cluster Distance = 0.118 



Tim e,t du = 9.51 


Fig . 5 Simulation for </>o(s) = 0 q(s). Left: ^(q) evaluated at Roq E /C. Right: Simulations of 
0 for t < 150. The vertical dashed line is at t = t c i u . 


dashed line) correctly estimates the time to the formation of the top and bottom 
clusters. 

We see an interesting result in Figure|4jfor </>o(s) = 0.5xlj o -y/y/ 2 ] ( s )+l(i/v^ i] ( s )- 
From the plot of we can see that q ma x might not be unique and the first few 
local maximizers are q = 2.5133, 9.4248, 15.7080,..., and the corresponding inter¬ 
cluster distances are 2.5, 0.6667, 0.4000,.... We can see from the simulation that 
there are two noticeable inter-cluster distances: 2.5 and 0.6667. For Roq £ JC : 

0 < q < 100, Qmax = 53.4071 so that the necessary condition to have the consensus 
convergence 9 max > 27 t is satisfied. However, we do not have consensus conver¬ 
gence in this case because <j m ax > 27r is not a sufficient condition. We notice that 
although (/max might not be unique, t c i u = 11.5 still predicts the time to cluster 
formation because it is related to ^(g m ax) not g ma x and thus the non-uniqueness 
of (/max does not affect t c i u . 

In Figure [5] and [6] we see consensus convergence. From the plot of tp(q), we can 
see that q m ax might not be unique and the first few local maximizers are q = 2.5133, 


































































12 


M 8 ) — 0.1 x 1 [0,1/V2] (s) + 1 X lpyy^l] (s) 
N = 500, Ro = 1, qmax = 53.4 
Inter-Cluster Distance = 0.118 
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Fig. 6 Simulation for 0o(®) — 0 q (-0 ■ Simulations of H for t < 1000. 



^o(s) = (1 - s) 3 X 1 [ 0 , 1 ] (s) 

N = 500, Ro = 1 , Qmax = 5.03 


Inter-Cluster Distance = 1.25 



Tim e,t du = 93.3 


Fig . 7 Simulation for (j>o(s ) = </>q(s). Left: "0(9) evaluated at Roq G /C. Right: Simulations of 
0- The vertical dashed line is at t — t c i u . 


9.4248, 16.3363,..., and the corresponding inter-cluster distances are 2.5, 0.6667, 
0.3846,.... In this case, the only noticeable inter-cluster distance is 0.6667 and we 
do not observe the inter-cluster distance of 2.5, because il>(2.5 133) <C ' 0 ( 9 - 4248 ). 
For Roq £ JC, 0 < q < 100, (/max = 53.4071 so that the necessary condition to have 
the consensus convergence (/max > 27r is satisfied and indeed we see form Figure [ 6 ] 
that we have consensus convergence in this case. 

In Figure [7j we choose 0o(s) so that 0(s) has a unique local maximum and 
(/max = 5.0265. In this case, (/max < 2 - 7 T and therefore there is no consensus conver¬ 
gence. The inter-distance is 1.25 and L/1.25 = 8 which is exactly the number of 
clusters in this case. Again, t c i u predicts the time to cluster formation very well. 

Finally, Figure [ 8 ] considers 0>o(s) so that ip(s) has a unique local maximum, 
but with a larger exponent, (/ ma x = 8.1681 > 27r. The inter-cluster distance corre¬ 
sponding to (/max = 8.1681 is 0.769, which is consistent with the actual inter-cluster 
distance. The quantity T/0.769 = 13 gives a good approximation for the actual 
number of clusters, which is 11. As in all the previous cases, t c i u = 185 predicts 
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M s ) = (! - s) 6 X 1 [o,i] (s) 

TV = 500, Rq = 1, Qmax = 8.17 


Inter-Cluster Distance = 0.769 



Fig. 8 Simulation for </>o(s) = 0 q(.s). Left: ip(q) evaluated at Rqq S /C. Right: Simulations of 
JTtI . The vertical dashed line is at t = t c i u . 


the time to cluster formation well. Here, g ma x > 27r, so we could expect to ob¬ 
serve consensus convergence. However the inter-cluster distance 0.769 is such that 
0(0.769) ~ 10 4 , so we cannot see cluster evolution for the time horizon of the 
simulation. 


5 Stochastic consensus convergence: er > 0 

In this section, we consider the case that a > 0 in ([TJ. In other words, the system 
is stochastic and we are dealing with a nonlinear Fokker-Planck equation. 


5.1 Linear stability analysis 


As in the deterministic case, we linearize the Fokker-Planck equation © with 
<j > 0 by assuming that p(t,x) = po + Pi(t,x) = 1/L + pi(t,x). Substituting 
p(t, x) = p 0 + pi ( t , x) into © and assuming that pi is a small perturbation of p so 
that the O(pf) term is negligible, we find that pi(t,x) satisfies: 


= po I ^(t,x-y)y(P{\y\)dy+^-^-(t,x). 

In the Fourier domain: 


dpi 

8t 


(t,k) 


ip 0 k j e lkv y(j>(\y\)dy 


2 i 2 1 

a k 
2 


pi(t,fc), 


(18) 


(19) 


which gives the growth rates of the modes: 


ipok J e lkv ycj>{\y\)dy - 



sin(ky)ycti(\y\)dy 


( 20 ) 


7 k = Re 


2 
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We can rewrite the growth rate 7 *. = poRo'J’aikRo), where 

i><r(q) = 2q J <po(s)s sin(qs)ds — ^ r3 ■ ( 21 ) 

The optimal positive frequency is fc m a X that is the element of K, = {27rn/L, n £ N} 
that maximizes il>a{kRo), that is close to qm&x/Roi where 


<?max = argmax [^< 7 ( 9 )] ■ 

q>0 


( 22 ) 


There is a critical value a c of <7 such that the system has a completely different 
overall behavior for a < a c and for a > a c - We can view a as the magnitude of the 
noise energy or temperature of the system. There are two types of forces in the 
system |lj): the attractive interaction —0^(27 — Xj) and the random force ad,W l (t). 
If a < a c , then the attractive interaction dominates the random force and thus the 
system is a perturbation of the deterministic system. If <7 > cr c , then the random 
force dominates, the attractive interaction is negligible, and therefore the overall 
system behaves like a system of N independent random processes. 

The above observations can be articulated mathematically. Let 


o> = max 
q> 0 


= max 

q>0 


4po Rfi 


1 fpo-Ro 


/ (po(s)s sin(qs)ds 

Jo 

f (/>o(s)s 2 Sln ^ qS ^ ds =Ap 0 Rl f s 2 (f>o{s)ds. 
Jo I s J Jo 


If <t < <7 C , then from (211 we find that 


maxi I>a{q) > 0, 

q> 0 


(23) 


and p{t, fc m ax) has positive growth rate 7max = Po-RoV’a-(^max-Ro)- The linear system 
( |19| ) is unstable, which is analogous to the deterministic case. 

If a > a c , then all of p(t,k ) have negative growth rates. In other words, the 
constant density is linearly stable and therefore the overall system is stable, since 
this is what linear stability implies in this case. 

We note that the same technique for computing a c , with linear stability analysis 
for different noisy opinion models, is also used in EMESES]. 


5.2 Fluctuation theory 

Since our goal is to analyze the behavior of clusters, we suppose from now on that 
a < a c - 

The fluctuation analysis of the stochastic model is similar to that of the de¬ 
terministic case. If xi(0),... ,®jv(0) are independent, uniform random variables in 
[0, L\, then 


( N 

pi (■t = 0, dx) := VN (p N (t = 0, dx) - po{dx)j = Vn f <5^ (0 )(rfm) 
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converges in distribution as IV — ► oo to the measure pi(t = 0, dx), whose frequency 
components 

= 0 ,k)= Ihn^ VN g e~ ik ^ j 

are independent complex circular Gaussian random variables, with mean zero and 
variance 1 for k £ Xi\{0}: 

E [pi(t = 0, k)] = 0, E (t = 0, k)pi {t = 0, £/)] = S kk ,, k,k'eK.\{ 0}, 


pi (t = 0, — k) — pi (t = 0, k), while pi (t = 0, k = 0) = 0. 
For any T < oo, the measure-valued process 


pi (t, dx) := y/N (p N (t, dx ) — po(dx)^) , f 6 [0, T] 


(24) 


converges in distribution as N —> oo to a measure-valued process pi(t,dx) whose 
density pi(f, x) satisfies a stochastic PDE (see [T]): 


dpi{t,x) = 


Po J d ^{t,x-y)y(f)[\y\)dy+ 


dt + adW(t,x) (25) 


with the random initial condition described above. Here W{t,x) is a space-time 
Gaussian random noise with mean zero and covariance 


Cov (J W(s,x)fi(x)dx, J W(t,x)f 2 (x)di 


iin{s, t} 


[ fi{x)f 2 (x)di 
Jo 


( 28 ) 

for any test functions /i (x) and /2 (x). The Fourier transform of W(t, x) is W ( t , k) = 
W(t,x)e~ lkx dx for k £ JC. From (26), we see that {W(t, k),k £ K,,k > 0} are 
independent, complex-valued Brownian motions with the variance: 

Cov (w(t, k), W(t, k)^j = J (— ik)e~ lkx {—ik)e~ lkx dx = tk 2 , (27) 

and W(t,—k) = W(t,k). Taking the Fourier transform on (25), for each k £ 1C, 
pi(t, fc) = f k pi(t,x)e~ lkx dx is a complex-valued Ornstein-Uhlenbeck (OU) pro- 


dpi(t,k) = 


rL 

ipok / e~ lkv ycf>(\y\)dy — 
Jo 


2 j 2 

a k 


PI (t, k)dt + °^d{W {k) ( t) + iW (k) (t)). 
v 2 


(28) 

Here {W (k \t),k £ AC, k > 0} and {W {k \t),k £ IC,k > 0} are independent real 
Brownian motions, W {0) (t) = W {0) (t) = 0, W ( ~ k \t) = -W {k) {t), and W { ~ k) (t) = 
W^ k \t). The equation (28) is solvable and we have, for any k £ A7\{0}, 

pi(t,k) = e“ fct pi(0,fc) + °JL [ e aUt ~ s) d{W ik \s)+iW {k \s)), (29) 

v 2 Jo 

a k = ip 0 k j e~ lky y<j>(\y\)dy - 
Jo 


2 
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In particular, pi(t,k = 0) = 0. Because {pi(0,k),W ( ' k \t),W < ' k \t),k £ K.,k > 0} 
are independent, {pi(t,k),k £ IC,k> 0} are independent OU processes with mean 
zero and variance 


E Jpi (i, k)p\ (t, fc)J = e 2 ' 1kt + a 2 k 2 J e 2lk< ' t S ^ds = e 2lkt 


1 + 


a 2 k 2 
2 7fe 


(30) 


for k £ 1C, k > 0, where 'yj- is the real part of a *.. In addition, because a_y = ay. 
and W { ~ k \t) = -W^ k \t), we have pi{t,—k) = p\{t,k) for k £ /C\{0}. Finally, 
pi(i,0) = 0. Therefore, 


E [p 1 (t,x)pi(t,x')\ = E 

= E E 


ikx 


E P 1 fa k ) ~r E p 1 (*>' 


Akx 


^ikx g— ikx' 

pi(t, k )~[^pi (*> _fe ) L 


_ _j_ g 27 k t e ik{x-x) 

r 2 


fc^O 


1 + 


27 fc 


As t increases, the spectrum of pi(i,a;) becomes concentrated around the opti¬ 
mal wavenumber fc m ax- In addition, we note that k 2 / 7 *. is bounded and the term 
a 2 k 2 / 7 *. is negligible if a is sufficiently small. We can assume a is small because we 
need a < a c for cluster formation. If a 2 k 2 /'yy is negligible and L 2 O 4 | 7 ^, ax |t, we 
can expand 7 /. = 7 max+ | 7 max(fc — &max ) 2 for k around fc max , and use a continuum 
approximation for the discrete sum as we do in the deterministic case: 


E [pi(t, x)p\(t, x)] 


(le 2 —‘ cos(fc max (a: 



( 1 

\ >/ ^iTmaxI^ 



If L 2 < 4|7max|^? then the continuum approximation is not valid and in this case 





cos 


(^max(^ X )) 


1 + 


2 » 2 

& ^max 
27max 


Because a < o c we have that 7 max > 0 and then the linear system ( |25| ) is 
unstable and therefore the central limit theorem breaks down when pi (t, x ) /y/N 
is no longer smaller than po = 1/L. More precisely, the time t c i u for the onset of 
clustering is when E[pi(i clu ,a;) 2 ] ~ NL~ 2 , which is approximately 


^clu — 


27, 


In A ~ -— „ , . -rlniV 

max ^PoRotycT (^max) 


when N 1. 


5.3 Consensus convergence 


We assume a < a c so that there are unstable modes for the linearized evolution, 
which means that there is clustering. The number of and the distance between 
clusters can be estimated with q ma x. We find that the first term of the right side 
of ( 21 ) is bounded while the second term of (211 is quadratic with negative leading 
coefficient. Therefore, increasing er tends to reduce g ma x, that is to say, to increase 
the mean distance between clusters. 
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Let us consider the case that q m ax < 2n. From the analysis of the deterministic 
case, the system initially has no consensus convergence and there are several clus¬ 
ters. After clustering, the clusters do not interact with each other, but their centers 
move like independent Brownian motions. When two clusters come close to each 
other, within a distance Rq, they interact and merge. Therefore, we will eventually 
have consensus convergence, because two Brownian motions always collide in R. 
This can be extended to the multi-dimensional case, but then the conclusion can 
be different: in high dimension two Brownian motions may not collide. However, 
with periodic boundary conditions, two Brownian motions will always come close 
to each other, within a distance Rq. 

When (j max < 27 t and a is small then there are several clusters, after the cluster 
formation time. The fraction of agents in a cluster is the agent density times the 
inter-distance of the clusters: 

2-7r.Ro 2-7T 

trie = po - = Pot -• 

Q'max ^max 

Then the j-th cluster consists of about Nm c agents. We assume that a is small 
enough so that cr 2 fc max 2irpo<j>o(0)Ro. By using the fact that the agents in a 
cluster stay close to each other, we can replace <j>(xi — Xj ) by </>(0), and the agents 
in the j-th cluster have the approximate dynamics: 

dx:p = — ( x^ — x\^)dt + adW^’ l \t). 

l=i 

The center X^\t) = satisfies: 

X u \t) = X (i) (0) + , a W u \t), (31) 

y/Nmc 

where {W^\t)} are independent standard Brownian motions. 

When N is large, the empirical density $ O) (dx) is approximately 

a Gaussian density 


p^\t,dx) = 


exp - 


{x-X^{t)f 


a (j) = 


/m c (j> o(0) 


For this argument to be valid, we must have that cr^\ the width of p^\ is much 
smaller than Rq , which is equivalent to our assumption <r 2 fc ma x 27rpo<( ) o(0)f?o. 

This cluster dynamics is valid as long as the centers {-X'^(t)} stay away from 
each other by a distance larger than Rq . The clusters move, according to indepen- 
dent Brownian motions with quadratic variation Nm . When two clusters come 
close to each other within a distance Rq, they merge. Indeed, once the two centers 
are within distance Rq, they obey the following differential equations to leading 
order in N: 


dX ^ = -m c (X (fe) (t) - X m (t))<i>(X ik \t) - X m (t)), 
dXi ^ ^ = -m c {X {l) (t)-X {k \t))<j>(X (k) {t)-X {l \t)), 
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which shows that the inter-cluster distance converges exponentially fast to zero 
and the center converges to the average of the centers just before collision. The 
number of agents or mass of the new cluster is the sum of the masses of the two 
clusters, the inverse square width of its empirical density is the sum of the inverse 
squares of the two widths, its center is at the weighted average (weighted by the 
masses) of the two centers just before collision and it moves as a Brownian motion 
whose diffusion constant is defined in terms of its new mass. Then the cluster 
centers move according to Brownian motions until two of them get within the 
distance Ro from each other and a new merge event occurs. This eventually forms 
a Markovian dynamics described in the next section. 


5.4 Markovian dynamics of the clusters 

After the initial clusters are formed, we can use an iterative argument to mathe¬ 
matically describe how all the opinions converge eventually. In the initial configu¬ 
ration, at time ro, there are M(ro) = Lk max /2n clusters with centers X^\tq) = 
j2n /Imax (up to a global shift), widths ct^(to) = a / ^/m c 4>o{ti), and masses 
m ^(r 0 ) = rn c for j = 1,... ,M (r 0 ). 

For t > t n— i, there are M(r n _ i) clusters moving as 

X (j \t ) = X U \T n -i) + a U \r n -(w^it) - kF (i) (r„_i)j 


until the stopping time 

T n = inf jf > r n - 1 : |X^(t) — X^ l \t)\ = Ro, for some k ^ l\ . 

Then the two colliding clusters (with indices k and l ) merge with the new center 

y, -k raW(T ra _i)fy( fc) (Tn) +m^(r ra _ i)N W (t~) 
m( fe )(r„_i) + m( i )(r n _i) 


the new mass 

m(r n ) = m (fc) (r„_i) + m (l) (T n -i) 


and the new width 


d(Tn) 


a 

v / m(r„)()>o(0) 


The clusters are relabeled to take into account this merging so that there are 
M(rn) = M(r n -i)—l clusters. The above process is repeated until n = (Lk m ax/27r) — 
1, when we have only one cluster, and hence consensus convergence. 

Note that the time scale at which collisions and merges occur is of the order 
of N, as the Brownian motions are scaled by 1/VN. 
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</>o(s) = l[o,i] (s),£T c = 0.36515,(7 = 0.1 
N = 500, Rq = 1, g ma x = 2.51 
Inter-Cluster Distance = 2.5 



Time, t du = 17.6 


Fig. 9 Simulations for a = 0.1. Left: 'ipa(q) evaluated at Roq G /C. Right: Simulations of 
The vertical dashed line is at t = t c [ u . 



5.5 Numerical simulations 


We use the explicit Euler scheme to simulate the stochastic opinion dynamics 0 
when a > 0: 

1 N 

Xi +1 =Xi - — -Xj\)(xi -x 1 j)At + aAW™ +1 , cf>(s) = cj> 0 (s/R 0 ), (33) 

i=i 

where {Vh™ +1 } are independent Gaussian random variables with mean zero and 
variance At. 


Our analysis is on the torus [0, L ], but we simulate (331 on [0, L] with reflecting 
boundary conditions. As we will see, the simulation results agree with the analysis 
under periodic assumption. Because here we focus on the effects of the randomness, 
for simplicity we will work only on the case that </>o(s) = 1 [o,i]( s )- 

We compute the key quantity g ma x by exploring all possible q in [0,100]: 


9max = argmax 

R 0 qelC,o<q<ioo 


' l 


2 q / (j>o(s)ssm(qs)ds — 


2 2 
o- q 


2pO-Rrj 


We see form the plots of ip<r(s) that the randomness reduces the possibility of the 
non-uniqueness of ax because it adds a negative quadratic term in ipcr{s). With 
randomness, all of our test cases have a clear, unique q max . 

The parameters we use for the simulation are At = 0.1, L = 10, Rq = 1 and 


N = 500. For each cr, we also plot the function ipa{s) in (21); the stars in the 
plots are the values of evaluated at Roq £ K, and the lines are the continuum 
approximation. 

We first to test for the effect of a c , the critical value for a, which makes the 
system stable or unstable. In our setting, a c = 4poAu fo s 2 4>o(s)ds = 0.365 and we 
simulate (33) for a = 0.1,0.2,0.365,0.5 that are values below, equal to and above 


(j c , respectively. 

From Figure [9j we see that i/>tr(q) decreases quadratically and has the unique 
maximum at q = 2.5133. However, maxjj oge jy ^(g) is still positive so the linearized 
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</>o(s) = l[o,i] (s),o- c = 0.36515,(7 = 0.2 
N = 500, R 0 = 1, gmax = 1.88 


Inter-Cluster Distance = 3.33 



0 50 100 150 

Time, t du = 33.9 



10 Simulations for a = 0.2. Left: evaluated at Rqq £ K,. Right: Simulations of 



system (191 is still unstable. Therefore, the overall system behavior is similar to 
the deterministic case and can be viewed as a perturbed non-random opinion 
dynamics. 

We increase a by setting am 0.2 and the result is in Figure |10| We see that 
as a increases, the random noise starts to affect the overall system, and the width 
and the inter-cluster distances become larger so we observe fewer clusters. Since 
max/j og6 x; ipcr{q) is positive, we still observe cluster formation. 

if a > (7 0 , il>a(q) < 0 for all q > 0 and 


We note that in Figure 


11 


and 


12 


Vv(0) = 0. In other words, the linearized system is stable and thus the full system 
is stable. In this case, we do not see cluster formation and the system behaves like 
an ^independent agent system. 

We see from the simulations that to model opinion dynamic with consensus 
convergence it is appropriate to assume that a < a c - Therefore we will assume 
that a = 0.1 in our simulations of the stochastic system. 

We revisit Figure [9] to check our analysis. First of all, g ma x = 2.5133 is clearly 
a unique maximizer and the corresponding inter-cluster distance is 2.51, which 
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0o («) = 1[0,1] (s), CT C = 0.36515,(7 = 0.5 



Time 



Simulations for a = 0.5. Left: ?/> CT (g) evaluated at Roq G 1C. Right: Simulations of 


agrees with the numerical inter-cluster distance we see in Figure [9] In addition, 
L/2.51 = 3.9841 also predicts well the actual number of the clusters, 4. Finally, the 
blue dashed line t = t c i u indicates the time to the cluster formation, even though 
Figure [9] is just one realization. 

We test t c i u and the width of clusters in a more statistical way by examining 
1000 realizations. If 0o(s) = l[o,i](s) and <r = 0.1, then we can expect that we will 
have 4 clusters at T = 150 in most of the realizations. For each realization, we 
numerically compute the widths of the clusters, (t), j = 1,..., 4 by using the 


empirical standard deviations of {x\ 3 \t)}j =1 (see (32)): 


(d«(t)) 2 = 


NU) - 1 


N u> 2 JV 0) 

Y ( 4 ^ w _ * 0) co) . x (3 \t) = y ^(t) ( 34 ) 


where for each j = 1,...,4, {abelong to the j-th cluster and is the 


number of agents in the j-th clusters. Of course, a 1 ' 3 ' 1 ( t ) in (34) is just one realization 
and so we compute a < - 3 ' 1 ( t ) for 1000 realizations and consider the average. 


13 


in different colors. 


The averages a^ 3 \t) are shown in the left part of Figure 
First, we can see that t c i u , as expected, is the halfway from the time to maximum 
with to the time to the minimum width. Second, from (32), the width of each 
cluster is analytically (t* j) = a/^/m c (j> o(0) = 0.1/V0.25 x 1 = 0.2, which agrees 
with the numerical values (t ) when t is large. 

We also analyze the behavior of the centers of the clusters. The centers {X ( T (t)}4 
of the clusters in Figure[9]are plotted in Figure [l4| From the previous analysis ( |31| |, 
the centers of the clusters are independent Brownian motions aW^ 3 \t)/s/Nm c . For 
one realization, the opinions {®j(t)}^Y_ 1 will not be evenly distributed in the clus¬ 
ters. For example, the actual numbers {N^}jL i of agents of the clusters in Figure 
[ 9 ] are plotted in Figure |l3| So for one realization, X^{t) is a Brownian motion 
aW^ 3 \t)/VW(j ). On the right part of Figure [m} we compare the quadratic varia¬ 
tions of X^ 3 \t) and aW^ ( t )/%/TVW) for 75 < t < 150 (after the time to the cluster 
formation.) Indeed, from the figure, we can see that their quadratic variations are 
very similar and that means X- 3 > (l) are very close to Brownian motions. 
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<t>0 (s) = 1[0,1] (s ) 5 cr = 0.1 
N = 500, Rq = 1 



12 3 4 

Clusters 


Fig. 13 Left: Average of 1000 realizations of j = 1,... ,4 for a = 0.1. The vertical 

dashed line is at t = t c i u . Right: the number of agents in the jf-th clusters in Figure [ 9 ] 



Time Time 


Fig. 14 Left: Centers of the clusters, X^(t), in Figured] in different colors. Right: Quadratic 
variations of the cluster centers from t = 75 to t = 150. Solid lines: Numerical quadratic varia¬ 
tions of xW(t). Dashed lines: Quadratic variations of ( t )/VJVW), where (fyO) (t)}4_i 

are independent standard Brownian motions and {-/V^}4 =1 are the numbers of agents in the 
clusters. 


6 Long time behavior of simulations 

We have also simulated numerically the long time behavior of the system defined on 
the full real line R, especially the behavior after the onset of consensus convergence. 
As we discuss in the previous section, when there is randomness the center of 
the unique cluster behaves like a diffusion process aW(t)/\/N, where W{t ) is a 
standard Brownian motion. In Figure [15] and [l6j we observe that the centers indeed 
behave like Brownian motions. The dashed lines are the parabolas with equation 
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<M S ) = 1 [ 0 , 1 ] ( s ); a = 0.1 <j>Q (s) = 1[0,1] (s), 0 - = 0.1 

N = 100, -R 0 = 1 N = 200,^o = l 



Fig. 15 Long time behavior of the opinions for a = 0.1 and for N = 100,200. The blue 
dashed curve is the equation x = d=2cr^/ t/N. a is small and the overall behavior is like a single 
Brownian motion. 


<A>0) = 1 [o,i]( s )i ,j = 0- 2 



M s ) = l[o,i] («), cr = 0.2 
N = 200, Ro = 1 

30 



Fig. 16 Long time behavior of the opinions for cr = 0.2 and for N = 100, 200. The blue 
dashed curve is the equation x = ±2cry/t/N. a is small and the overall behavior is like a single 
Brownian motion. 


x = ±2a^t/N so that for any fixed f, the centers are within the parabolas with 
95% probability. 

However, when a is sufficiently large, the long time behavior is different. On the 
right part of Figure[l7] when a = 0.4 > <r c , the system behaves like ^-independent 
diffusions. A more interesting case is when a = 0.3 < cr c on the left part of Figure 


for t > 5.5 x 10, all Xi(t) spread out from the unique cluster and the system 
becomes an independent agent evolution. A detailed mathematical analysis using 
large deviations theory for such a phenomenon is being considered at present. 


17 In this case, for 0 < t < 5.5 x 10 5 there is still consensus convergence, 
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M s ) = 1 [ 0 , 1 ] (s), cr = 0.3 Ms) = l[ 0 ,l]( s )> a = 0.4 

N = 100,R 0 = 1 N == 100, Rq = 1 



Time * 10 ‘ Time 


Fig. 17 Long time behavior of the opinions for a = 0.3, 0.4 and for TV = 100. The blue dashed 
curve is the equation x = ±2u\/ t. For a = 0.3, there is a single cluster for t < 5.5 X 10 5 . 
However, for t > 5.5 X 10 5 , Xi(t) disintegrate and the system behaves like an TV-independent 
Brownian motions. For a = 0.4, the random perturbations are large enough so that the system 
is an TV-independent Brownian motions at the beginning. 


7 Conclusion 


We have analyzed a stochastic, continuous time opinion dynamics model and we 
have carried out extensive numerical simulations. We use the mean-field theory 
and obtain a nonlinear Fokker-Planck equation as the number of opinions tends 
to infinity. Then we use a linear stability analysis to estimate the critical value 
of the noise strength so as to have cluster formation, estimate the number of 
clusters and the time to cluster formation. These quantities are closely related to 


the frequency that maximizes the growth rate of the linearized modes (201. After 


the initial cluster formation we expect, and numerically confirm, that the centers 
of the clusters behave like Brownian motions before further consolidation. Finally, 
the long time behavior of the system is explored numerically and we observe that 
after a unique cluster is formed, there is a small probability that the opinions will 
spread out from the unique cluster and the system will become an independent 
agent evolution. 
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